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ABSTRACT 

We consider the structure of marginally Toomre-stable starburst disks under the assumption that radiation pres¬ 
sure on dust grains provides the dominant vertical support against gravity. This assumption is particularly ap¬ 
propriate when the disk is optically thick to its own infrared radiation, as in the central regions of Ultralumi- 
nous Infrared Galaxies (ULIRGs). We argue that because the disk radiates at its Eddington limit (for dust), the 
“Schmidt-law” for star formation changes in the optically-thick limit, with the star formation rate per unit area 
scaling as E* oc Eg/ k, where Eg is the gas surface density and k is the mean opacity of the disk. Our calculations 
further show that optically thick starburst disks have a characteristic flux, star formation rate per unit area, and 
dust effective temperature of F ^ 10'^ Lq kpc“^, E* ~ 10^ Mq yr“' kpc“^, and Teff ^ 90 K, respectively. We 
compare our model predictions with observations of ULIRGs and find good agreement. 

We extend our model of starburst disks from many-hundred parsec scales to sub-parsec scales and address the 
problem of fueling active galactic nuclei (AGN). We assume that angular momentum transport proceeds via global 
torques (e.g., spiral waves, winds, or a central bar) rather than a local viscosity. We consistently account for the 
radial depletion of gas due to star formation and hnd a strong bifurcation between two classes of disk models: (1) 
solutions with a starburst on large scales that consumes all of the gas with little or no fueling of a central AGN and 
(2) models with an outer large-scale starburst accompanied by a more compact starburst on 1 - 10 pc scales and a 
bright central AGN. The luminosity of the latter models is in many cases dominated by the AGN, although these 
disk solutions exhibit a broad mid- to far-infrared peak from star formation. We show that the vertical thickness 
of the starburst disk on parsec scales can approach r, perhaps accounting for the nuclear obscuration in some 
Type 2 AGN. We also argue that the disk of young stars in the Galactic Center may be the remnant of such a 
compact nuclear starburst. 

Subject headings: galaxies:general — galaxies:formation — galaxies:starburst — accretion, accretion disks — 

Galaxy: center — quasars:general 

1. INTRODUCTION 

Star formation in galaxies is observed to be globally ineffi¬ 
cient; in a sample of local spiral galaxies and luminous star- 
bursts, Kennicutt (1998) showed that only a few percent of the 
gas is converted into stars each dynamical time. This ineffi¬ 
ciency may result from “feedback”: the energy and momentum 
injected into the interstellar medium (ISM) by star formation 
can in turn regulate the star formation rate in a galaxy. Models 
for feedback in the ISM generally invoke energy and momen¬ 
tum injection by supernovae and stellar winds (e.g., McKee & 

Ostriker 1977; Silk 1997; Efstathiou 2000). However, the mo¬ 
mentum supplied to the ISM by the radiation from massive stars 
is comparable to that supplied by supernovae and stellar winds. 

The UV radiation from massive stars is absorbed and scattered 
by dust grains, which reprocess the UV emission into the IR. 

Because the dust grains are hydrodynamic ally coupled to the 
gas, radiation pressure on dust can help stabilize the gas against 
its own self-gravity and may therefore be an important feedback 
mechanism. 

When the ISM of a galaxy is optically thick to the re-radiated 
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IR emission, radiative diffusion ensures that all of the momen¬ 
tum from the photons produced by star formation is efficiently 
coupled to the gas. We show that this limit is applicable on 
scales of hundreds of parsecs in luminous gas-rich starbursts, 
including Ultraluminous Infrared Galaxies (ULIRGs), the most 
luminous and dust enshrouded starbursts known (e.g., Genzel 
& Cesarsky 2000). Indeed, Scoville (2003) has already pointed 
out that radiation pressure on dust could plausibly be the domi¬ 
nant source of support against gravity in ULIRGs. We quantify 
this hypothesis by developing models of radiation pressure sup¬ 
ported starburst disks. 

Although radiation-pressure supported disks have not been 
extensively considered in the galactic context, they have been 
well studied in models of black hole accretion. Most notably 
for our present purposes, the outer parts of disks around active 
galactic nuclei (AGN) are expected to be dominated by radia¬ 
tion pressure on dust (e.g., Sirko & Goodman 2003). In the con¬ 
text of this paper, the distinction between the galactic disk and 
the AGN disk becomes somewhat unclear: if luminous AGN 
are fueled by gas from the cold ISM of their host galaxy, there 
must be a continuous transition from the star-forming “galac- 
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tic” disk to the central black hole’s “accretion disk.” The nature 
of this transition, and indeed whether it occurs at all, remains 
uncertain. The problem is that the outer parts of AGN disks 
are strongly self-gravitating with a Toomre stability parame¬ 
ter Q ^ 1 (Kolykhalov & Sunyaev 1980; Shlosman & Begel- 
man 1989; Shlosman, Begelman, & Frank 1990; Kumar 1999; 
Goodman 2003; Levin 2005; Tan & Blackman 2005). It is dif¬ 
ficult to see how the disk avoids fragmenting almost entirely 
into stars. One possibility is that in the dense gas-rich nuclear 
regions of galaxies, angular momentum transport proceeds via 
global torques (e.g., bars, spiral waves, large-scale magnetic 
stresses, etc.), rather than a local viscosity (Shlosman, Frank, 
& Begelman 1990; Mihos & Hernquist 1996; Goodman 2003). 
In this case, gas may inflow sufficiently rapidly to avoid turning 
entirely into stars. In order to explore the fueling of a central 
AGN, we extend our galactic-scale disk models to smaller radii 
using a parameterization of the transport of angular momentum 
by global torques. This provides a framework for understanding 
the relationship between AGN activity and nuclear starbursts. 

The plan for this paper is as follows. In §2, we develop a sim¬ 
ple model of self-regulated star formation in galactic disks. We 
distinguish between disks that are optically thick/thin to their 
own IR radiation and we present models appropriate to each 
of these limits. Throughout §2, and in particular in §2.5, we 
compare our theoretical models with observations of starburst 
galaxies. In §3, we extend the results of §2 to include angu¬ 
lar momentum transport and we address the problem of AGN 
fueling. Finally, §4 provides a discussion and summary of our 
conclusions, and highlights some of the predictions of our mod¬ 
els. 


2. A ONE-ZONE DISK MODEL 

In this section, we construct simple dynamical models for the 
structure of star-forming disk galaxies. We consider two limits: 
(1) the “optically-thin” limit, in which the disk is optically-thick 
to the UV radiation produced by massive stars, but optically- 
thin to the re-radiated IR radiation and (2) the “optically-thick” 
limit, in which the disk is optically-thick to the re-radiated IR. 

We begin by describing the properties of our models com¬ 
mon to both limits. We make the clear oversimplification that 
the disk is a single phase medium (see, however, the brief dis¬ 
cussions in §2.1 and Appendix B); our model thus describes 
only the average properties of the disk. We assume that the 
disk is in radial centrifugal balance with a rotation rate 12 = 12 k, 
where 12 k = v^cr/r is the Keplerian angular frequency in an 
isothermal potential with velocity dispersion a. The total dy¬ 
namical mass at radius r is given by Miot(r) = 2cP-r jG and the 
associated surface density is 

cP' 

Stot = ~ 0.6 ct^oo ''kpc g cm■^ (1) 


that this model captures much of the physics of interest on large 
scales. In §3 we relax this assumption and show that a constant 
gas fraction is equivalent to a constant accretion rate. 

For a thin disk in a spherical potential, the equation for ver¬ 
tical hydrostatic equilibrium, dp/dz = -pVP'z, can be approxi¬ 
mated as 

( 2 ) 


where h is the pressure scale height. This implies h « Cs/12 
where cj = p /p 'k the sound speed. Throughout this paper we 
include turbulent pressure in our definition of the sound speed 
(in addition to radiation pressure and gas pressure). A term 
l-KGTigph should be added to the right-hand-side of equation 
(2) to account for the self-gravity of the disk. This leads to a 
multiplicative correction to equation (2) of (1 -l-2^/^/0, which 
we neglect for simplicity. 

We assume that star formation in the disk is governed by 
local gravitational instability as described by Toomre’s Q- 
parameter (Toomre 1964). In particular, we argue that the disk 
self-regulates such that 

rtoCo 12^ 

is maintained close to unity. In equation (3), = 4f2^ -I- 

c/fl^/c/lnr is the epicyclic frequency. The hypothesis of 
marginal Toomre stability has been discussed extensively in the 
literature (e.g., Paczynski 1978; Gammie 2001; Sirko & Good¬ 
man 2003; Levin 2005) and is based on the idea that if I 
then the disk will cool rapidly and form stars, while if Q<^\ 
then the star formation will be so efficient that the disk will heat 
up to 2 ~ 1. There is evidence for Q ~ 1 in the Milky Way (e.g., 
Binney & Tremaine 1987), local spiral galaxies (e.g., Martin & 
Kennicutt 2001), and starbursts such as ULIRGs (e.g., Downes 
& Solomon 1998). 


From equation (3) it follows that the density distribution of 
the gas is determined solely by the local Keplerian frequency: 


f2^ 


TxGQr^ 


170cr|ooi-kpcG ‘cm ^ 


(4) 


We retain the Q dependence here and below for completeness, 
but our assumption is that 2 ~ 1. From the definition of and 
equations (1) and (4), it then follows that 


(/r/r) = /,e/23/2 (5) 


and that 

{c,la) = fgQl2. ( 6 ) 


For constant fg and Q, hjr and Cj/ct are independent of radius. 
With cr = 200 km s“', a turbulent sound speed of 10 km s“' cor¬ 
responds to a gas fraction of , = /g/0.1. The scale height of 
the disk for this gas fraction is then 


where (T 200 = cr/200kms“^ and the radial scale is in units of 
kpc. For simplicity, we assume that the underlying potential 
is spherical, as would be provided by a stellar bulge or the 
galaxy’s dark matter halo.^ We further assume that the gas mass 
is a constant fraction fg = = S^/Etot of the total dy¬ 

namical mass, where Eg is the gas surface density. Although 
the assumption of a constant fg is clearly simplistic, we believe 


/*~354o,i'‘kpcPc. (7) 

In our model, the star formation rate per unit area E* adjusts 
to maintain 2 ~ 1. We also parameterize star formation as oc¬ 
curring on a fraction rj of the local dynamical timescale (e.g., 
Elmegreen 1997): 


E* = Egf2p. (8) 

^ For a Navarro, Frenk, & White (1997) dark matter profile, the dark matter mass can be significant even on the small scales of interest: M(r) ^ this 

result is relatively independent of halo mass so long as r where r? is the scale radius of the NFW potential. 
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Only values of r; < 1 are physical. For rj greater than unity the 
disk cannot dynamically adjust to maintain Q ~ 1. Although, 
S* is the fundamental derived quantity, 77 is a useful alternative 
parameter that characterizes the global star formation efficiency 
in the disk. Observationally, p is typically ~ 0.02 in normal spi¬ 
ral and starburst galaxies (Kennicutt 1998). 

2.1. The Optically-Thin Limit 

We consider galactic disks to be “optically-thin” when the 
vertical optical depth to IR photons is < 1, but the optical 
depth to UV photons is greater than unity. This requires 
2 /kuv ~ 10 ^ g cm where kuv ~ 10 ^ cm^ g ' is a character¬ 
istic UV opacity. 

In our models, feedback from star formation provides the 
pressure support to maintain Q ~ 1. Sources of this pressure 
support include radiation pressure on dust grains, supernovae, 
stellar winds, and expanding HII regions. In this section we 
consider a simple model that relates the pressure in the ISM to 
the star formation rate E* in the optically thin limit (eq. [15] 
below). To motivate this model, and to connect with classic 
treatments of the ISM (e.g., McKee &. Ostriker 1977), we esti¬ 
mate the pressure in both the “cold” {pc) and “hot” (p/,) phases 
of the ISM, focusing on the contributions from radiation pres¬ 
sure and supernovae. The pressure in the cold ISM can thus be 
written as: 

= ~Prp + Psn, (9) 


Supernovae also generate a shocked hot ISM, whose pres¬ 
sure ph can be estimated using the model of McKee & Os¬ 
triker (1977). Individual supernovae initially expand adiabat- 
ically to a cooling radius Rc ~ 20Ef^n~^ ‘^ pc, where a cool 
dense shell forms (Chevalier 1974; Cioffi et al. 1985; Thorn¬ 
ton et al. 1998). The subsequent remnant evolution is approx¬ 
imately momentum-conserving. In the absence of interaction 
with neighboring supernova remnants, the remnant expands to 
a maximum radius Rp ~ PC, where the pres¬ 

sure of the supernova remnant is equal to the ambient pressure 
(P) and P -12 = P/10“'^ ergs cm“^. In fact, the remnant typically 
does not reach Rp before encountering a neighboring supernova 
remnant. That is, the overlap radius for supernova remnants 
{Ro) is Rc < Ro < Rp- In this model, the average pressure of 
the hot ISM is not known a priori, but can be estimated using 
the pressure of a supernova remnant at Rg (this assumes that 
there is no significant energy loss after the supernova remnants 
overlap and mix). Using results from Chevalier (1974), McKee 
& Ostriker (1977) showed that the pressure at overlap can be 
expressed in terms of the pressure at Rc as* 

Ph~Pc^'c, (12) 

where = SVcL is the the number of supernovae that occur 
in the volume Vc ~ ( 4 / 3 ) 7 rP*, C is the time for the remnant to 
reach Rc, and S is the supernova rate per unit volume. The time 
tc has been estimated by Cox (1972) and Cioffi et al. (1985): 

fc ~ 4 X yr for solar metallicity. Using Pc = E jVc 

we find that 


where q should be interpreted as a turbulent velocity and pc as 
a turbulent pressure. 

Because we assume that the optical depth to UV photons is 
greater than unity, the radiation pressure on cold gas can be re¬ 
lated to the star formation rate by 

Prp ~ eS*c ( 10 ) 

where e is the efficiency with which star formation converts 
mass into radiation (e ~ 10“* for a Salpeter IMF from 1 - 
IOOMq). Equation (10) assumes that there is no significant 
cancellation of oppositely directed momentum. This becomes 
an increasingly better assumption as the disk becomes opti¬ 
cally thick to the IR (as in starburst galaxies). In addition 
to radiation pressure, supernova explosions deposit momen¬ 
tum into the ISM via swept-up shells of cold gas. From 
Thornton et al. (1998) we estimate that each supernova has an 
asymptotic momentum of 3 x g cm s“', where 

£51 = £/ 10 ^^ ergs is the initial energy of the supernova and 
n\=nl\ cm“* is the density of the ambient medium. The contri¬ 
bution to the ISM pressure from supernovae may be estimated 
by taking psn = An/(27rr^), where Am is the total momentum 
injection rate. Alternatively, one can estimate psn by balanc¬ 
ing the net volumetric energy injection rate by supernovae, 
~ AnCj/(27rA/i), with energy losses due to turbulent de¬ 
cay. Assuming that turbulence decays on a crossing time so 
that ~ p(?slh (e.g.. Stone et al. 1998), we see that 
implies that Am/(27rr^) = = psn- Thus, the contribution to 

the pressure of the cold ISM from supernovae is 

Psn ~ 1.5 X ~ 5n7‘/%‘i^/^Vrp (H) 

Equations (10) and (11) show that radiation pressure and super¬ 
novae contribute comparably to the pressure of the cold ISM. 


Ph ~ 10 ergs cm (13) 

where we have scaled the number density to a value comparable 
to that of the hot ISM (n _2 = n/10“^ cm“*), and the supernova 
rate to the Galactic value: 5 -i 3 = 5/10“'* pc“* yr“*. Equation 
(13) is in reasonable agreement with observations and more de¬ 
tailed models of the local pressure of the ISM (McKee & Os¬ 
triker 1977; Boulares (& Cox 1990). 


In order to compare p\^ and prp from equation (10) directly, 
we note that 5 ~ 10 “*E*/ {2K) with E* in Mq yr“' per unit area. 
The ratio of these two components of the pressure can then be 
written as 


Prp 


3 / 7 -' 7,-4/V 
jni00”-2 ^51 > 


(14) 


where /,ioo = /z/lOO pc. Because the total volume occupied 
by remnants with R < Rc decreases as the density of the ISM 
increases, the contribution from supernovae to the total pres¬ 
sure decreases with increasing density. Equations (9)-(14) show 
that, to order of magnitude, pc = Prp + Psn ~ Ph for conditions 
appropriate to the Galaxy (as is observed). 


In the luminous starbursts we focus on in this paper, the den¬ 
sity of the ISM is much larger than in the local ISM. Eor exam¬ 
ple, in the inner few hundred parsecs of ULIRGs the average 
gas density reaches 10 *- 10 ^ cm“*, comparable to the density 
of a local molecular cloud (e.g., Downes & Solomon 1998; see 
also eq. [4]). Several lines of evidence suggest that the cold 
molecular gas may fill a significant fraction of the volume in 
ULIRGs, unlike in the local ISM (Downes, Solomon, & Rad¬ 
ford 1993; Solomon et al. 1997). The high fraction and high 
luminosity of radio supernovae in Arp 220 is also consistent 
with an environment much denser than the ISM of normal spiral 


More carefully, should be raised to the 10/11-th power in equation (12), which we approximate as unity for simplicity. 
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galaxies (Smith et al. 1998). Taking n ~ 10^ cm“^ as a charac¬ 
teristic value, we find that Rc is just ~ 1 pc, p ^/~ lO'^/zj'o, 
and that the total asymptotic thermal energy of a supernova 
remnant is ~ 4 x 10'** ergs <C E (Thornton et al. 1998). This ar¬ 
gues against a dynamically-dominant, volume-filling hot ISM. 
Even in the limit of strong radiative losses, however, super¬ 
novae are still important for generating the random motion of 
cold gas (Psn ~ ^’rp > Ph\ eq. [11]). 

In the simple estimates above, all of the contributions to the 
pressure of the ISM scale roughly linearly with the star for¬ 
mation rate. Moreover, the ratio p^/p^„ is of order unity, and 
while Ph/Prp is of order unity in normal spiral galaxies it may be 
significantly smaller in the dense nuclei of the most luminous 
starbursts. For this reason we choose to express the effective 
pressure of the ISM in the “optically-thin” limit as 

P~Pe=Prp[l+(Psn/Prp)] = eCS*C, (15) 

where the last equality defines the parameter In what follows 
we retain the dependence of our results on but scale to ^ ~ 1 
for the reasons given above. 


Using equations (4), ( 6 ), and (15), it is straightforward to 
solve for the physical parameters of our disk model. The star 
formation rate per unit area required to support the disk with 
2 ~ 1 is given by 

_ fya 

Scaling equation (16) for typical fg and a we find that 

~ 354 \,cr^oo'‘kpce 3 'r‘G Mo yr”' kpc “2 

where €3 = e/10“* and fg^^ = fg/Q.5 is appropriate for gas-rich 
starbursts. The star formation rate can also be expressed in 
terms of the efficiency p; 


(16) 

(17) 
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0-l/go,5'^200e3'^ * 


(18) 


The second equality in equation (18) is meant to show explicitly 
that one may write rj in terms of Cj or fga (eq. [6]). Using equa¬ 
tion (16), the flux and luminosity of the disk viewed face-on are 
given by 


F = ~ 5 X Lq kpc'^ (19) 

and 

^^ ^"" io'VL'"2ooer‘ Lq. (20) 

Up to logarithmic corrections, equation (20) implies that for 
constant fg all radii contribute equally to the total luminosity. 


The effective temperature is defined by the relation ctsb = 
eE*c^. In the optically-thin limit, the observed dust temperature 
Tdust is related to Teff by ~ 7’gff/(2rv), where ry is the ver¬ 
tical optical depth to IR radiation. This relation between Tjust 
and Teff assumes that the sources of UV radiation are uniformly 
distributed vertically throughout the disk. With equation (16), 

Tiu. = (^^372!^) - 604 />^/o'(MCrkpc)-'/^G^/^K, 

( 21 ) 

where ki = k /\gem ^ is a representative value for the IR opac¬ 
ity of the disk (Fig. 1) and we have assumed that an individual 
dust grain radiates as a blackbody. 


There are several interesting properties of the disk model pre¬ 
sented in this section. First, the efficiency of star formation 
required to maintain G ~ 1 is p ~ 0.02 for a canonical tur¬ 
bulent velocity of Cs ~ lOkms”* (equivalent in our model to 
having fg ^ 0.1 for a ~ 200 km s as in the Galaxy). This 
value for rj is in good agreement with observations compiled 
by Kennicutt (1998). Second, the first two equalities in equa¬ 
tion (16) yield s* cx This Schmidt-like star formation law 
is somewhat steeper than the S* cx favored by Kennicutt 

(1998), but comparable to the scaling S* cx obtained by 
Gao & Solomon (2004) using a sample that includes more lumi¬ 
nous starburst galaxies. Given the simplicity of the model pre¬ 
sented here, this agreement is satisfactory. Third, for fg = 0.1 
and cr = 200 km s“* the dust temperature (eq. [21]), turbulent 
velocity, pressure, flux, luminosity, and scale height (eq. [7]; 
compare with Fig. 9.25 from Binney & Merrifield 1998) are all 
in fair agreement with observations of the Milky Way. Thus de¬ 
spite its simplicity, the model presented in this section provides 
a useful characterization of galactic-scale star formation sup¬ 
ported by the turbulent pressure produced by supernovae and 
the radiation from massive stars. 

2.2. The Optically-Thick Limit 

The nuclei of gas-rich starbursts are optically thick to their 
own infrared radiation. Radiative diffusion then ensures that ra¬ 
diation pressure provides the dominant vertical support against 
gravity. In this section we describe disk models appropriate to 
this limit. The vertical optical depth of the disk is given by 
Ty = EgK/2, where k is the Rosseland mean opacity to dust. 
Evaluating this expression yields 

W = ~ 0-15CT2oo/«o.5'-kpc«i- (22) 

The radius at which ry = 1 is then 

^rv=l = - 150KiCrfoo4„, pC (23) 

and thus for the largest most gas-rich starbursts (cr ~ 300 km 
s“* and fg ~ 1) the inner ~ 700 pc are optically thick. 

In the optically thick limit, the effective temperature is given 
by 

cr SB Trff = — eS*c^, (24) 

where the factor of 1 /2 arises because both the top and bottom 
surfaces of the disk radiate. The midplane temperature is re¬ 
lated to the effective temperature by « (3/4)ry7’gff, where 

the opacity k{T,p) in ry should be evaluated using the central 
temperature and mass density of the disk. The temperature de¬ 
pendence of the opacity is important for this problem and is 
discussed in the next section. For now we simply normalize k 
to 1 cm^ g“'. 

For Ty > 1, the radiation pressure is given by 

4crsB_4 CTsb 1 ^ 

Frad = —T = -ryTgff = -TveS*c. (25) 

3c c 2 

Comparing equations (15) and (25) shows that radiation pres¬ 
sure exceeds the turbulent pressure due to supernovae by a fac¬ 
tor of ~ Ty (assuming ^ ~ 1 for the reasons given in §2.1). Ra¬ 
diation pressure support will thus dominate the vertical support 
of compact optically-thick starbursts. The exact surface density 
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(or Tv) at which the transition to radiation pressure support oc¬ 
curs is somewhat uncertain and requires a better understanding 
of the pressure support provided by supernovae, stellar winds, 
etc.. We will explore this in more detail in future work. 


With equation (25) it is again straightforward to solve for 
the disk’s physical parameters, T, T^ff, E*, etc. in terms of our 
model variables: a, fg, e, and the radius in the disk r. The 
midplane temperature of the disk is 


3c6 


and the effective temperature is 

1/4 

t S III \ 

Teff = 


1/4 


fgo cQ \ 
r V2KcrsB / 


W41a2oo4/^kpfe'/^K. (26) 


(27) 


For our fiducial numbers, Teff is somewhat larger than T, im¬ 
plying that the vertical optical depth is less than unity. Our as¬ 
sumption that the disk is optically thick in the far IR is therefore 
only marginally applicable on kpc scales. It is an increasingly 
better assumption on smaller scales (see eqs. [22] and [23]). 


Therefore, the “Schmidt-law” for star formation changes in the 
dense optically-thick inner regions of starburst galaxies. 

Equation (31) can be rewritten as 
L _ F _ nQGcfg _ 2 ttGc h ^^3 fg,,, Q Lq 

M Etot ^/2K nr k\ M@ 

The third equality above is the classical Eddington limit, mod¬ 
ified by the factor {h/r) ex fg, and with the Rosseland mean 
opacity taking the place of the usual electron scattering opac¬ 
ity.® This way of expressing the flux highlights the fact that 
each disk annulus radiates at its local Eddington limit. The 
value derived here for the luminosity per unit mass from the 
disk is similar to that estimated by Scoville (2003), who ar¬ 
gues that this limit is observed in both young star clusters such 
as M51 and ULIRGs such as Arp 220 (see also Scoville et 
al. 2001). Note that for a given gas fraction, the mass-to-light 
ratio is proportional to the dust opacity and is therefore metal- 
licity dependent. 


2.3. Opacity Dependence 


The total star formation rate per unit area required to support 
the disk with radiation pressure is 

~ 400/g„5cr|oo2'‘kpc«T‘e3‘ M© yr“' kpc“l 

T 

(28) 

Integrating, we derive the total star formation rate: 

2^/^7rf Oa^r 

M* =- ^ -3OOO/g„5(J2oo?-kpc0e3'KT' M© yr“^ (29) 

Note that if the disk is optically-thick (2 x tv ~ 1), the observed 
flux is only from half of the disk and therefore the observation- 
ally inferred star formation rate would be one-half of the true 
total star formation rate (given in eqs. [28] & [29]). 

The star formation efficiency rj is 

nGQ r _i _i _i 

p =-~ 1 rkpcCr20oe3 Q- (30) 

/CCC (7 

Eor r < Rtv=i we see that p is < 1. It follows that the disk can 
adjust to maintain 2 ~ 1 on sub-kpc scales. In addition, we find 
that the ratio of the cooling timescale to the orbital timescale is 
much less than unity, so the disk should self-regulate to main¬ 
tain 2 ~ 1 (see Appendix B). 

Einally, the surface brightness of the disk viewed face-on is 

F = ~ 3 X 10'V,„.5<^20o<e«T' L© kpc-2 (31) 

v2k r 

and the total luminosity for a single side of the disk is 

^ = '^'^fgQcro- ^ 2 X lO^V«o.5O'2Oo'‘kpc0K7' L©. (32) 

K 

Equations (27)-(32) show that the observable properties of 
the disk depend sensitively on the magnitude of the opacity. In 
particular, the star formation rate per unit area is proportional 
to E^/tv ~ Sg/tv (eq. [28]). Therefore, in regions of the disk 
where the opacity is low, there must be more star formation to 
maintain Q ^ 1. Conversely, where the optical depth is high, 
less star formation is required. The functional dependence of 
E* in the optically-thick limit should compared with the equa¬ 
tion (16), which shows that in the optically-thin limit E* oc E^. 


Eigure 1 shows the Rosseland mean opacity as a function 
of temperature for several densities using the publicly avail¬ 
able opacities of Semenov et al. (2003) and Bell & Lin (1994). 
There are two important features. Eirst, for temperatures T < 
100-200 K, the opacity is essentially independent of density 
and can be approximated by k = kqT^, with kq ~ 2.4 x 10“^ 
cm^ g“' K“^. The scaling of n with follows from the fact 
that the dust absorption cross section scales as with 5 —> 2 
in the Rayleigh limit (Pollack et al. 1985). The normalization 
Ko is somewhat uncertain and depends on grain physics and the 
dust-to-gas ratio; the latter may vary systematically as a func¬ 
tion of radius and metallicity in starburst disks. In fact, in our 
own galaxy there is evidence that the dust-to-gas ratio increases 
within the central few kpc (Sodroski et al. 1997). In what fol¬ 
lows we set k_3,6 = ko/2.4 X 10“® cm^ g“' and retain the 
scaling with kq. The second important feature of Eigure 1 is 
the dramatic decrease in the opacity for 10^ K < T < 10® K, the 
“opacity gap.” Here, the temperature is larger than the subli¬ 
mation temperature of dust, but smaller than the temperature at 
which hydrogen is significantly ionized. 


Equation (26) shows that, even for the largest galaxies with 
the highest gas fractions, the temperature at ^0.1-1 kpc is 
< 100 K, and so the opacity on large scales can be approxi¬ 
mated by K = kqT^. We may then eliminate the opacity de¬ 
pendence from the disk properties derived above. Remarkably, 
because T ex and E* cx E^/ k, with «; oc we find that the 
star formation rate per unit area, the effective temperature, and 
flux are all independent of virtually all model parameters: 


E* = 


/ 2®/VGasBe \ 
V 3€^kIc^ ) 


1/2 

~ M©yr-‘kpc-2 


(34) 


F = 


2^l^'KGa5BQc \ 
3^0 ) 


1/2 

-- L©kpc“^ 


(35) 


and 


Teff = 


f2^I^TTGQc\ 

V 3ko(Jsb J 


1/8 


88K_3;6e‘/*^ K. 


(36) 


' The factor 27r rather than the more standard 47r in equation (50) is because the luminosity in equation (32) is only from 1/2 of the disk. 
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In particular, note that neither S*, F, nor depend on r, a, 
or fg. Our model thus predicts that starburst disks have roughly 
constant flux and effective temperature over a range of radii. 

The constancy of these disk observables follows from three 
ingredients: (1) the disk is supported by radiation pressure and 
Tv > 1, (2) the disk self-regulates with 2 ~ 1, and (3) k (xT^ 
at low T. Above T ^ 100-200 K, K ceases to increase mono- 
tonically with T (see Fig. 1) and equations (34)-(36) no longer 
hold. Because T oc (eq. [26]), the temperature exceeds 
~ 200 K at a radius /?200 ~ ^OctIoo/^os^iooG^^ pc- This radius 
should be compared with the radius at which Ty = 1. Using the 
K (xT^ scaling, we find that 

R^=i = a^fg (^ 2 wZ2(7tG)^ ) ^ 

(37) 

(compare with eq. [23]). Therefore, we expect S*, Teff, and F 
to be roughly constant in the radial range R 20 Q < Rrv=i, ~ 
hundreds of parsecs for fiducial parameters. 

For completeness we note that the scaling for the dust tem¬ 
perature in the optically-thin limit can also be re-written using 
k<xT^ (eq. [21]): 

Tdust ~ (38) 

The weak scaling with model parameters in equation (38) im¬ 
plies that the dust temperature should not vary significantly 
from system to system, as appears to be observed (e.g., Yun 
& Carilli 2002). 

2.4. Combining the Optically-Thin & Optically-Thick Limits 

The optically-thin and optically-thick limits can be combined 
by expressing the pressure as 

p = eE*cQTv-l-C^ (39) 

using the full temperature-dependent opacity curve (Fig. 1). 
Recall that in the limit ry ^ 1 equation (39) describes true ra¬ 
diation pressure support while in the limit Ty ^ 1 it describes 
turbulent support with contributions from supernovae and ra¬ 
diation pressure. Because the opacity depends on temperature, 
we must connect the central temperature with the effective tem¬ 
perature in order to solve equation (39). By interpolating be¬ 
tween the optically-thin and optically-thick regimes we obtain 
(see also Sirko & Goodman 2003) 

In the optically-thick limit, equations (39) and (40) combine 
to yield (4(TsB/3c)r‘* ~ p, whereas in the optically-thin limit 
(4a-sB/3c)2Tv7’‘* ~ p. In solving equations (39) and (40), we 
find multiple solutions because of the complicated temperature 
dependence of the opacity. In Appendix B, we address the ther¬ 
mal and viscous stability of these solutions. We argue that there 
is a single stable physical low-temperature solution and focus 
on this solution throughout the paper. 

Figure 2 shows the numerically calculated structure of our 
disk models for cr = 200 km s“* and with fg=0.03 and 1. There 
are three regimes to notice in Figure 2. First, at large radii the 
disk is optically thin. In this region S*cxr 2 and ry ~ const (see 
§2.1). At smaller radii the disk becomes optically thick. There 


is then a range of radii (50 < r < 300 pc for fg = l) where 
E* and Teff are roughly constant, in good agreement with the 
estimates in §2.3. Note that where the disk is optically thin 
^dust ~ 7].|f/(2'fv)- At very small radii ~ 1 - 10 pc the opacity 
decreases dramatically when dust grains sublimate (the “opac¬ 
ity gap”; Fig. 1). In this region the disk becomes optically thin, 
and the star formation rate required to maintain Q ~ 1 increases 
significantly (see also Sirko & Goodman 2003). We present a 
detailed discussion of this part of the disk in §3 and Appendix 
A, but note here that it is unphysical to assume that the gas frac¬ 
tion is constant throughout the region where the star formation 
rate increases so markedly. 

2.5. Application to ULIRGs 

To focus the discussion, we emphasize the application of our 
optically thick disk models to Arp 220, a prototypical ULIRG. 
Arp 220 consists of two merging nuclei separated by about 350 
pc (Graham et al. 1990). The total FIR luminosity of the sys¬ 
tem is ~ 10 '^ Lq. The 2- 10 keV X-ray luminosity is only 
~ 3 X 10^ Lq, however, and the column density of X-ray ab¬ 
sorbing material must exceed 10^^ cm“^ if an obscured AGN 
is to contribute significantly to the bolometric luminosity (Iwa- 
sawa et al. 2001). Thus, there is little evidence for an energet¬ 
ically important AGN. The detection of numerous radio super¬ 
novae also supports a starburst origin for most of the radiation 
from Arp 220 (Smith et al. 1998). 

The stellar velocity dispersion of Arp 220 is ~ 165 km s“' 
(Genzel et al. 2001) and it has an extended CO disk with a 
scale of ~ 500 pc (Downes & Solomon 1998). Perhaps ~ 1/2 
of the luminosity, however, appears to come from two counter¬ 
rotating nuclear disks, each of which is ~ 100 pc in extent (e.g., 
Sakamoto et al. 1999; Soifer et al. 1999; Downes & Solomon 
1998). Radio observations also show that nearly all of the radio 
flux — presumably associated with supernovae and star for¬ 
mation — originates within ~ 50- 100 pc of the double nu¬ 
clei (Condon et al. 1991). The total mass in each of the com¬ 
pact nuclear disks is ~ 2 x IO^Mq, with a large gas fraction 
of fg ~ 0.5 (Downes & Solomon 1998). The nuclear region 
of Arp 220 is optically thick to at least 25/rm and the inferred 
blackbody temperature at this wavelength is > 85 K (Soifer et 
al. 1999). The estimated gas mass and radius of the disk imply 
a surface density of Eg « 10 g cm“^, suggesting that the nuclear 
region is probably optically thick even in the FIR (as implied by 
eq. [23]). 

In the optically-thick limit, with k oc T^, our model predicts a 
radiative flux of ~ lO^^LQkpc"^ (eq. [35]), in good agreement 
with the flux inferred for the compact nuclei in Arp 220 (see 
also Fig. 3 discussed below). Our model also predicts an effec¬ 
tive temperature of ~ 88 K for this very compact emission (eq. 
[36]), in agreement with that estimated by Soifer et al. (1999). 
Lastly, we note that for a gas fraction of fg ~ 0.5 and a velocity 
dispersion of cr ~ 200 km s“\ our model predicts a sound speed 
(i.e., turbulent velocity or line-width) of Cs ~ 50 km s“*, similar 
to that observed in CO by Downes & Solomon (1998). 

In their study of CO emission from ULIRGs, Downes & 
Solomon (1998) identified compact nuclear starbursts in Mrk 
273 and Arp 193 with properties similar to those discussed 
above for Arp 220. More directly, Condon et al. (1991) im¬ 
aged the radio emission in the 40 brightest galaxies in the IRAS 
Bright Galaxy Sample. They resolved the emission in nearly all 
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of the sources (36/40) and found sizes ~ 100 pc. They further 
showed that the radio emission from ULIRGs — correcting for 
free-free absorption in some cases — places them on the FIR- 
radio correlation of starburst galaxies. This is consistent with 
starbursts contributing significantly to the bolometric power of 
these sources. In the absence of direct FIR imaging of ULIRGs, 
the radio sizes from Condon et al. (1991) are currently the best 
estimate of the size of the nuclear starbursts in these systems. 
Figure 3 shows a histogram of the number of sources at a given 
flux inferred using the FIR luminosity and the size of the re¬ 
solved radio source from Condon et al. (1991). The histogram 
shows a peak centered on a flux ~ lO'^ Lq kpc“^. This charac¬ 
teristic ‘observed’ flux is consistent with the predictions of our 
optically thick disk models (eq. [35]).^® 

An alternative way to present this data is shown in Figure 
4, where we plot the inferred flux as a function of the size of 
the resolved radio source. Superimposed on the data we plot 
the predictions of our disk models for cr = 200 km s~^ and sev¬ 
eral fg (solid lines), and for a = 300km s~* and fg = l (dashed 
line). The latter is a plausible upper limit to the emission ex¬ 
pected in our starburst models. This Figure demonstrates the 
excellent agreement between our models and the observations 
for fg ~ 0.3- 1. The data are also consistent with the increase 
in flux we predict for more compact starbursts. Taken together, 
the results of Figures 3 and 4 provide strong support for our 
interpretation of ULIRGs as Eddington-limited starbursts. 

It is also worth comparing these results to an empiri¬ 
cal surface brightness limit for starbursts found by Meurer 
et al. (1997), who argued that starburst fluxes sat¬ 
isfy ^ 2 X 10^' L© kpc“^, which corresponds to E* < 
13 € 3 ' Mq yr“' kpc“^. Our calculations, and the observations 
summarized above, suggest that Meurer et al.’s limit does not 
represent a fundamental limit to the surface brightness of star- 
burst galaxies. 

The characteristic flux ~ 10^^ Lq kpc“^ found in Figure 3 
is equivalent to a blackbody temperature of ~ 90 K (eq. [36]). 
This is noticeably larger than the typical color temperature of 
~ 60 K inferred from the FIR spectra of ULIRGs. Another way 
to state this result is that using the observed FIR spectra and lu¬ 
minosity, the blackbody size of the FIR emitting region is typi¬ 
cally larger (by a factor of few) than the radio sizes observed by 
Condon et al. (1991). This is likely because the compact nuclei 
of many ULIRGs are optically thick even at ~ 30/im, so that 
radiative diffusion degrades the ^ 90 K emission and ensures 
that the FIR size can be larger than the true size of the nuclear 
starburst. This interpretation requires sufficient obscuring gas 
at large radii, but also that star formation in this gas does not 
dominate the bolometric power of the source (or else the ra¬ 
dio source would be more extended than is observed). In our 
models we find that if the gas fraction in the nuclear region in¬ 
creases at small radii, as Downes & Solomon (1998) infer for 
several systems, then most of the luminosity is produced near 
the radius where Ty ~ 1, rather than in the extended optically 
thin portion of the disk at larger radii. 

To conclude this section, we compare our model predic¬ 
tions with the observations of high-redshift ULIRGs (“submm 
sources”) presented in Tacconi et al. (2005). Tacconi et al. de¬ 
scribe CO observations of 4 high-z systems with the PdBI. They 

Note that our conclusions are not significantly changed if only (say) ~ 1/2 
some models of Arp 220 suggest; e.g., Soifer et al. 1999). 


find an average velocity dispersion of cr ~ 290 ±35 km s“\ 
gas fraction of fg ~ 0.4 ± 0.2, and half-power radius for CO 
of Rqo — 1.6 ± 0.3 kpc. The global properties of the disks in 
the high-z ULIRGs appear to be scaled up versions of their lo¬ 
cal counterparts, with luminosities ~ 10 x higher and CO-disks 
~ 3 X larger, implying comparable fluxes and E*. For the fg 
and a inferred from the observations, we predict that the ra¬ 
dius at which the disk becomes optically-thick to its own IR 
radiation is Rtv=i — 400-500 pc. Our disk calculations yield a 
total luminosity of 1 - 2 x 10'3 Lq and a flux that varies from 
6 X 10 '^ Lq kpc“^ at 1.6 kpc to lO'^ Lq kpc“^ at 200 pc. We 
also predict a dust temperature of ~ 57 K at ~ 1.6 kpc, a sound 
speed of Cj ~ 60 km s“', and a gas surface density of Eg ~ 0.8 
g cm“^. These predictions are in excellent agreement with the 
observations of Tacconi et al.: L ~ 10^^ Lq, F ~ IO'^Lq kpc“^, 
Cj ~ 95 ± 30 km s“^ and Eg ~ 0.6 ± 0.1 g cm“^. As discussed 
extensively above, in local ULIRGs the nuclear starbursts ap¬ 
pear significantly concentrated with respect to the large-scale 
CO disks. It remains to be seen if the same is true in their high- 
redshift counterparts (Chapman et al. 2004 find radio sizes of 
~ 3 kpc in several systems, suggesting that star formation may 
be more extended in some high redshift ULIRGs). 


3. DISK MODELS WITH ACCRETION 


In the previous section we considered the properties of large- 
scale starburst disks with constant gas fraction. In this section 
we address the problem of AGN fueling by connecting our kpc- 
scale starburst disks with AGN disks on sub-parsec scales. This 
requires a consistent treatment of the gas fraction, which must 
evolve with radius as a result of star formation. 


On large scales, models with constant gas fraction are equiv¬ 
alent to models with a constant mass accretion rate. To see 
this, consider a Shakura-Sunyaev accretion disk with viscosity 
= acsh. The mass accretion rate in such a disk is given by 


M = 2'Ki'Tig 


c/lnU 


c/lnU 

d\nr 

GQ 

c/lnr 


(41) 


Equating the scale height determined by equation (41) with that 
defined in equation (5) yields a one-to-one correspondence be¬ 
tween the gas fraction and the mass accretion rate: 


fg ~ 


( — 

c/lnU 

\aQ^ 

c/lnr 


MG \ 


1/3 


(42) 


Since U oc in an isothermal potential, equation (42) shows 
that a constant gas fraction implies a constant accretion rate. On 
small scales, where the black hole dominates the gravitational 
potential, fg oc for constant M. 

In reality, the accretion rate (or fg) is not constant. As gas 
accretes from larger to smaller radii and some of the gas is con¬ 
verted into stars to maintain 2 ~ 1, the accretion rate decreases 
monotonically. At the outer radius of the disk Rout, we assume 
that gas is supplied at a rate Mout- At any radius r < Rout the 
accretion rate M{r) is 

M(r)=Moui-J^ 27rr'S*c/r'. (43) 

of the bolometric power is produced by the compact nuclei resolved in the radio (as 



If the star formation rate required to maintain Q ~ 1 is large 
enough, all of the gas goes into stars and the accretion rate be¬ 
comes vanishingly small. To quantify this it is useful to con¬ 
sider two timescales characterizing the disk: (1) the advec- 
tion timescale Tadv = and (2) the star formation timescale 
T* = 1 /(rjil). At any radius in the disk, there is a critical accre¬ 
tion rate Me for which Tadv = Taking the optically-thick limit 
and assuming n = kqT^ (§2.3) yields 


Mc = r^ 


/ 2‘^/^7r^GQcrsB \ 


1/2 


(44) 

where we have scaled the radius to 200 pc in anticipation of nu¬ 
merical calculations described below. If the gas accretion rate 
in the disk satisfies M > Me, then Tadv < and the gas is able 
to accrete inwards to smaller radii. By contrast if M < Me, then 
Adv > 2* and most of the gas is converted into stars at radius r 
without flowing significantly inwards. We note that gas supply 
rates exceeding Me seem plausible in view of the fact that many 
ULIRGs have star formation rates greater than Me. 


Equation (44) has several important properties. First, Me 
is seemingly independent of the viscosity. This results from 
the fact that we have assumed k oc T^, which is valid only 
at low temperatures and thus large radii. For general k or in 
the optically-thin limit. Me depends explicitly on the viscosity 
(see Appendix D). Second, Me is an increasing function of ra¬ 
dius. Thus if gas is supplied at the outer radius Rout at a rate 
Mout > Me, significant gas accretion can continue to smaller 
radii. However, on scales smaller than ~ 10 pc, the central disk 
temperature is high enough that dust sublimates, n decreases 
dramatically, and the star formation rate required to maintain 
\ increases (Fig. 2; see also Sirko & Goodman 2003). The 
large star formation rate at small radii makes it difficult to fuel a 
central accretion disk at a rate sufficient to explain bright AGN. 
This competition between star formation (t*) and inflow (Tadv), 
particularly throughout the opacity gap on 0.1 - 10 pc scales, 
determines the rate at which a central AGN is fueled. 


large radii (eq. [44]), the fate of gas at small radii in the opac¬ 
ity gap is a very strong function of the rate of angular momen¬ 
tum transport. The efficiency of angular momentum transport is 
important because, for fixed M, higher viscosity implies lower 
surface density and thus self-gravity is comparatively less prob¬ 
lematic. Indeed, it is well-known that local angular momentum 
transport — such as is produced by the magnetorotational in¬ 
stability (Balbus 2003) — is incapable of supplying sufficient 
gas to a central black hole to account for luminous AGN (e.g., 
Shlosman & Begelman 1989; Shlosman et al. 1990; Goodman 
2003). 


One possible solution to this problem is that angular mo¬ 
mentum transport proceeds by global torques such as would be 
provided by stellar bars, spiral waves, or large-scale magnetic 
stresses (Shlosman et al. 1990; Goodman 2003)." In this sec¬ 
tion, we use a phenomenological prescription to describe this 
process: we assume that the radial transport of gas by a global 
torque allows the radial velocity to approach a constant fraction 
m of the local sound speed (Goodman 2003). In this case. 


M = 2TTrTigVr 




(45) 


and the relationship between accretion rate and gas fraction is 
given by 


p 3/2 mg 
\ Qm (Or)3 J 


(46) 


instead of equation (42). Our hope is that, much as the Shakura- 
Sunyaev prescription provides a useful zeroth-order model for 
local angular momentum transport in disks, the above model 
captures some of the essential physics of disks in which angu¬ 
lar momentum transport is dominated by global torques. With 
equation (45) to relate the gas surface density to the gas accre¬ 
tion rate, we solve the equations of Appendix C to determine 
the structure of the disk. 


3.2. Results 


3.1. A Model of AGN Fueling 

In Appendix C we collect the equations and parameters used 
in this section to derive the properties of starburst and AGN 
disks. As in §2, we assume that H = Hk, but here we ac¬ 
count for the gravitational potential of the black hole: = 

GMsh/+ 2a^I. The black hole mass is assumed to be given 
by the M^a-a relation: Mbh — 2 x lO^cr^oo ^0 (Tremaine 
et al. 2002; Ferrarese & Merritt 2000; Gebhardt et al. 2000). 
In calculating vertical hydrostatic equilibrium, we include gas 
pressure, which is important for r < 1 pc. The most important 
change relative to the model of §2 is the use of a consistent 
accretion rate that accounts for the loss of gas locally to star 
formation (eq. [43]). 

Although the critical accretion rate Me at which Tadv = 
does not depend that strongly on the viscosity in the disk at 


The dashed lines in Figure 5 show the mass accretion rate M 
as a function of radius for Mout = 80, 160, 220, 320, and 640 
M© yr“' in a model with a = 300 km s“' (Mbh — 10^ Mq), 
Rout = 200 pc, and with dynamical angular momentum trans¬ 
port specified by m = 0.2. The black hole dominates the gravi¬ 
tational potential for r < 50 pc. The solid lines show the local 
star formation rate M* defined by M* = Trr^S*." 

There are two types of solution represented in Figure 5. The 
two models with the smallest M^^t lose essentially all of the 
supplied gas at large radii to star formation. At r ~ 20 - 40 
pc in these two models the star formation rate becomes so low 
that M approaches a constant near ^0.1 M© yr“^ This oc¬ 
curs when the central temperature decreases sufficiently that 
gas pressure dominates. These models have Tadv > 0* at Rout 
and M^ut < MeiRoud- They are starburst-dominated and the star 
formation occurs predominantly at ^ Rout- 


' * An alternative possibility not considered here is that AGN are fueled by low angular momentum gas (Goodman 2003), including perhaps the hot ISM in clusters of 
galaxies (e.g., Nulsen & Fabian 2000). 

If fg exceeds unity at Rout, then we expect the disk to become globally self-gravitating. Through equation (46), this condition sets a “maximum” mass accretion 
rate, ~ 28 Mq yr“*, above which the approximations of this paper break down. For these reasons, we do not consider models with Mout > 640 

Mq yr“' for the parameters of Figure 5. Equation (5) also shows that as fg —> 1, h/r —> 1. For models in which h/r or fg exceeds unity, we expect m or a to 
increase because of the large gas fraction (as in the “bars within bars” model of Shlosman et al. 1989), a wind to be driven from the disk surface (thereby limiting 
Mout dynamically), or perhaps the gas to be pushed out radially (Rout increases) to maintain h/r ^ 1. 
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For the three models shown with Mout > 220 Mq yr ' the re¬ 
sults are qualitatively different. In these models, Mout is large 
enough that star formation persists for more than a decade in 
radius from /?out = 200 pc down to r ~ 1 - 10 pc. At these 
small scales, the temperature is sufficiently hot that dust sub¬ 
limates and the opacity decreases sharply (see Figure 1). Be¬ 
cause E* (X Ej/k in the optically-thick limit, the star formation 
rate must increase dramatically to maintain 2 ~ 1. The increase 
in M* as the opacity gap is encountered is typically an order of 
magnitude. Although these star formation rates are large, they 
are much smaller than those computed for the constant M (or 
constant fg) solutions in Figure 2 (§2). In the models with vari¬ 
able fg, an equilibrium between advection and star formation, 
expressed by the equality Tadv = rv, limits the star formation rate 
in the opacity gap (Appendix A). Each of these high Mout disk 
models has a gas accretion rate of ~ 4 M© yr“' at r ~ 0.01 pc, 
sufficient to power a luminous AGN. That these models have 
nearly identical accretion rates at small radii follows from the 
fact that Tadv = T* in the opacity gap. In Appendix A we present 
an analytical model for the structure of the disk in this region 
and provide simple formulae for the black hole accretion rate. 

As Figure 5 shows, there is a strong bifurcation between 
models that fuel a central AGN and those that do not. The crit¬ 
ical mass supply rate distinguishing these two classes of solu¬ 
tions is close to Me as estimated in equation (44) by equating 
Tadv and r*.'^ 

For comparison with these calculations assuming dynamical 
angular momentum transport, we have also calculated disk so¬ 
lutions employing a local a-viscosity (eq. 41). With a = 0.3, 
the disk structure exterior to the opacity gap is very similar to 
the models with dynamical angular momentum transport be¬ 
cause h/r is large on scales of /?out in the models with large 
Mout- Therefore, a local viscosity can transport gas from several 
hundred parsec scales down to ~ 10 pc scales in gas-rich star- 
bursts. However, we find that a local viscosity cannot transport 
gas through the opacity gap at a rate sufficient to fuel bright 
AGN; typical black hole accretion rates are several orders of 
magnitude smaller than those with dynamical angular momen¬ 
tum transport (see Appendix A). 

In the three models with Mout > 220 Mq yr“' > Me, we 
can distinguish between the outer “starburst disk”, where the 
heating of the disk is dominated by star formation (crsBTgg = 
(l/2)eE*c^), and the inner “AGN disk”, where the heating is 
dominated by accretion (crsBTgg = (3/87r)M(l - 
The transition between these two regimes occurs at /?agn ~ 0.5 
pc in the models presented in Figure 5. At this radius, accretion 
heating is sufficient to produce Q>\ and star formation ceases. 

By combining the equations describing the “starburst” disk 
at r > /?agn with those describing the “AGN disk” at r < /?agn 
we have computed the full radial structure of models from 
hundreds of parsecs down to the central black hole. We note 
again that there are multiple solutions at some radii and that 


we have selected what we believe is the physical disk solution 
as described in Appendix B. As an example of the full radial 
disk structure. Figure 6 shows the results for the model with 
Mout = 320 Mq yr“' (see also Fig. 5). In the middle left panel, 
note that gas pressure dominates radiation pressure over ~ 1 
decade in radius inside the opacity gap at r < 1 pc (Appendix 
A), but that radiation pressure dominates at all other radii. 

Assuming that every annulus in the disk radiates as a black- 
body and that the disk is viewed face-on, the spectral energy 
distribution of our disk models can be computed using 

lirhc^ /■*“' lirrclr 

^ jRin exp[/zc/AkB7;ff(r)]-r 

The resulting multi-color blackbody spectra for the models pre¬ 
sented in Figure 5 are shown in Figure 7, assuming an inner disk 
radius of Ri^ = 3(2GMbh/ c^)- The models with small mass sup¬ 
ply rate (Mout < Me) have the lowest bolometric luminosity and 
virtually all of their flux is generated by the starburst at Rqui- 
The spectra of these models peak at ^ 50 /rm, corresponding to 
an outer disk dust temperature of ~ 70 K. For higher Mout, the 
dust temperature at Rqui increases somewhat and the spectral 
peak moves to shorter wavelengths. 

Higher Mout models have similar FIR spectra. However, at 
~ 8 /im they exhibit an additional spectral peak. This peak is 
produced by the inner ring of star formation on 1 -10 pc scales 
shown in Figure 5, which is caused by the decrease in the opac¬ 
ity of the disk at the dust sublimation temperature, ^ 1000 
K. Although Tsub sets the central temperature of the disk where 
this sudden increase in E* occurs, the disk is typically optically 
thick at these radii and the effective temperature can be a fac¬ 
tor of ~ 3 - 4 smaller, accounting for the fact that the spectral 
peak associated with star formation in the opacity gap is at ~ 8 
pm. Together, the FIR peak from the outer starburst ring and 
the peak from the inner starburst ring create a broad IR peak 
spanning ~ 1.5 decades in wavelength. 

The distinguishing feature of the models with large Mout is 
the significant contribution of the AGN to the bolometric lumi¬ 
nosity: the UV emission in Figure 7 is comparable to or dom¬ 
inates the peaks at ~ 8 and ^ 50 pm. Figure 7 suggests that 
quasar accretion disks can plausibly be fed by radiation pres¬ 
sure supported starburst disks at larger radii. Indeed, the broad 
but sub-dominant IR peak and the magnitude of the AGN lumi¬ 
nosity in our models are reasonably consistent with composite 
quasar spectra (Elvis et al. 1994; Haas et al. 2003). 

It should be noted that the spectra shown in Eigure 7 are pho- 
tospheric and do not include AGN reprocessing. Because the 
disk may intercept some of the AGN emission, we expect that 
reprocessing could lead to an additional emission component at 
~ 1 - 2pm if dust in the surface layers of the disk is heated to 
the sublimation temperature. 


3.3. Vertical Structure & Nuclear Obscuration in AGN 


As an aside we note that all models presented in Figure 5 have small M* at small radii and one may worry about discreteness; that is, for very small M* only 
a few massive stars may be formed in a given annulus and because an individual massive star only supports a volume ~ around it the disk may not be capable 
of self-regulation. The “continuum” approximation employed throughout this work is only valid when r/h <C Nm*{r) ~ M*(r)rm*10“*'^^MQ, where is the 

number of massive stars at radius r, Tm* is the lifetime of a massive star, and we have assumed a Salpeter IMF and defined “massive stars” as those with mass greater 
than ~ 10 Mq. For Tm* in the range 10® — lO’ yr the continuum approximation is valid for the models displayed in Figure 5. For a discussion of when the continuum 
limit breaks down, see footnote 15. 

In these calculations, we have assumed that dynamical angular momentum transport acts throughout the disk. On sub-parsec scales near the black hole, perhaps 
at i?AGN where Q becomes greater than unity, we expect accretion to be driven instead by local viscosity. This will modify the density, scale height, etc., close to the 
black hole (< 0.5 pc in Figure 6), but will not change our conclusions about the need for angular momentum transport by global torques at larger radii. 
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At radii ~ 1 pc the models shown in Figure 6 have the inter¬ 
esting property that the central temperature of the disk is above 
the sublimation temperature of dust while the effective (surface) 
temperature of the disk is below the sublimation temperature. 
As a result, there is a large vertical gradient in the opacity of 
the disk. Because the disk is radiation pressure dominated, this 
large gradient in the opacity has important implications for the 
disk’s vertical structure. Here we argue that the photosphere 
of the disk can lie a distance hpk ~ r off of the midplane even 
though the midplane scale height is small, with h <^r (where 
h = Csj H is evaluated using the midplane properties of the disk). 
We then discuss the implications of this result. 


The equations governing the vertical structure of the disk in¬ 
clude those of vertical hydrostatic equilibrium and vertical en¬ 
ergy transport: 



(48) 


where p = pg^s + Prad and we assume that H is independent of 
the height z in the disk. Assuming that energy is transported 
diffusively by photons implies 


dT 3 kpF 

dz Ihr^CTSB 


(49) 


where F is the vertical flux of energy. If the disk is radiation 
pressure dominated, as our models are at most radii, then equa¬ 
tions (48) and (49) imply that the flux in the disk must be ev¬ 
erywhere equal to the Eddington flux: 


F = Fedd = 


n'^cz 

K 


(50) 


The presence of a large vertical opacity gradient suggests that 
convection might develop in the outer atmosphere of the disk. 
However, two arguments show that convection is unlikely to 
be important in the present context: (1) Our disk models have 
TthO 1, where t± is the cooling time of the disk (Appendix 
B). This implies that the diffusion time through the disk is much 
shorter than the characteristic timescale of buoyant convective 
motions. In the presence of such strong radiative diffusion, con¬ 
vective modes driven by radiation entropy gradients grow very 
slowly, if at all (e.g., Blaes & Socrates 2003). (2) Even if con¬ 
vection develops, it is unlikely to carry a significant fraction of 
the vertical energy flux. Using F « pV^ to estimate the convec¬ 
tive velocity Vc required to carry the energy flux we find very 
high Mach numbers: 


Yi 

Cs 


40M, 


1/3 



1/3 




-2 ^- 1/3 ^ 1/6 
F-15> 


(51) 


where p_i5 = p/10“^^ g cm“^ and we have scaled our estimate to 
parameters appropriate to parsec scale-disks (see Eig. 6). These 
considerations strongly suggest that radiative diffusion domi¬ 
nates the vertical transport of energy throughout the disk. 


Unfortunately, equations (48) and (49) are insufficient to 
fully specify the vertical structure of the disk. An equation for 
the vertical generation of energy, dF/dz, is also needed (e.g., 
Davis et al. 2005). The vertical heating profile is particularly 
uncertain in the present context where star formation, rather 
than accretion, dominates the heating of the disk. Given this 
significant uncertainty, we restrict ourselves to the following 
simple considerations about the disk’s vertical structure. 


We assume that the calculations described in the previous 
sections — that apply to the midplane of the disk — provide an 


adequate estimate of the star formation rate S* and, thereby, the 
flux F required to maintain Q ~ 1. This assumes that the atmo¬ 
sphere of the disk does not contribute much flux, a reasonable 
assumption. Because the disk is radiation pressure dominated 
near its photosphere, we can estimate the location of the photo¬ 
sphere using F « Fedd- This yields 


^ph ~ 


F Kph 



(52) 


where Kph is the opacity at the photosphere, evaluated using the 
effective temperature, and Kmid is the opacity at the midplane, 
evaluated using the central temperature. Equation (52) shows 
that if the opacity at the surface of the disk is much larger than 
the opacity near the midplane, as will inevitably occur in the 
vicinity of the sublimation temperature of dust, then hph ^ h 
and the midplane scale height does not provide a good estimate 
of the location of the photosphere of the disk. 


Eigure 8 shows our estimates of the midplane scale-height 
h/r and the photospheric scale-height hph/r for models with 
Mout = 320 and 640 M© yr“' (whose other properties have 
been described in Eigures 5-7). These estimates of the pho¬ 
tospheric scale-height show that hph can readily approach ^ r 
at ^ 0.1 - 10 pc, even though the midplane scale-height of the 
disk satisfies h r. Because the effective temperature of the 
disk is < Fsub, the opacity near the photosphere is k ~ 1 g cm“^ 
and the hydrogen column through the disk’s atmosphere can 
approach ~ 10^“* cm^. We suggest that this extended dusty at¬ 
mosphere may account for the presence of obscuring material 
inferred from observations of AGN (e.g., Antonucci 1993). In 
particular, although AGN disks are typically classified as being 
“thin” because h<^r, Eigure 8 shows that the photosphere may 
nonetheless puff up substantially and reach hph ~ r. Some of 
this material is likely to be unbound and the origin of a dusty 
outflow. 


The disk’s extended atmosphere is supported by radiation 
pressure from a parsec-scale nuclear starburst (in contrast to 
Pier & Krolik 1992, who argued that dusty gas could be sup¬ 
ported at /z ^ r by radiation pressure from the AGN). Estimating 
the luminosity of the disk at these radii, L = nr^F, we find 

^ 2rG/lfc /hph\ ^ Fgs / tves \ ^hph\ (53) 

Kph \ r J 4 V Kph y V r y ’ 

where Fes is the canonical Eddington luminosity defined using 
the electron-scattering opacity Kes. Taking Kph ~ IOkes, equa¬ 
tion (53) implies that in order to have hph ~ r, as suggested 
by the relative number of Type 1 & Type 2 Seyfert galaxies, 
the luminosity of the parsec-scale disk must be ^ 0.01 -0.1 Fes- 
Therefore, in order for the disk to have hph ~ r, we predict the 
presence of a very compact nuclear starburst with a luminosity 
similar to that of the AGN. 

Because the luminosity from star formation required to “puff 
up” the disk is substantial, it is unlikely that this model is rel¬ 
evant to low-luminosity Seyfert 2 galaxies such as NGC 4258 
(where there is also strong evidence that the nuclear obscuration 
is due to a warped disk; e.g., Eruscione et al. 2005). A compact 
starburst may, however, dynamically support obscuring mate¬ 
rial in more luminous AGN. In particular, there is evidence for 
a compact starburst in the prototypical Seyfert 2 Galaxy NGC 
1068. IR interferometry of 1068 has also directly resolved 
warm dusty gas on the parsec scales predicted by our model 
(Jaffe et al. 2004; Rouan et al. 2004). 
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4. DISCUSSION & CONCLUSIONS 

In Standard models of galactic-scale star formation, energy 
and momentum injected by supernovae and stellar winds are 
assumed to drive turbulent motions in the interstellar medium 
of a galaxy (e.g., Silk 1997). This turbulent pressure helps stave 
off the self-gravity of the disk and maintain marginal stability to 
gravitational perturbations (Toomre’s Q ~ 1). In this paper, we 
have focused instead on the role of radiation pressure on dust 
grains in regulating the structure and dynamics of star forma¬ 
tion in galaxies. To order of magnitude, the turbulent pressure 
from supernovae, stellar winds, and the radiation from massive 
stars are all comparably important when the galactic disk is op¬ 
tically thin to its own IR radiation, as in normal star-forming 
galaxies. By contrast, when the disk is sufficiently optically 
thick to the IR, radiation pressure provides the dominant ver¬ 
tical support against gravity. This condition is met in the in¬ 
ner few hundred parsecs of luminous gas-rich starbursts, most 
notably in Ultraluminous Infrared Galaxies (ULIRGs). In addi¬ 
tion, the outer parts of accretion disks around AGN are expected 
to be dominated by radiation pressure on dust. Understand¬ 
ing the dynamics of disks in the radiation pressure dominated 
limit is therefore crucial for understanding nuclear starbursts 
and AGN fueling. 

We have constructed simple quantitative models of self- 
regulated disks appropriate to both the optically-thin and 
optically-thick limits. From these models, we derive the star 
formation rate per unit area required to maintain 2 ~ 1. In 
the optically-thin limit, we find that S* ex (eq. [16]), a 
slightly steep version of the Schmidt law for star formation, 
whereas in the optically-thick limit, E* oc E^/ry (x E^/k (eq. 
[28]). Because radiation pressure dominates the vertical pres¬ 
sure support in the optically-thick limit, the star formation rate 
is sensitive to the mean opacity of the disk, and thus also to 
the metallicity. In fact, each annulus of the disk radiates at its 
local Eddington limit, defined using the Rosseland mean opac¬ 
ity K. This criterion can be written as (see also Scoville 2003) 
L/M~ 10Vg„.«T' Lq/Mq, where the gas fraction fg is propor¬ 
tional to the disk thickness h/r. 

Our prediction that the Schmidt law changes in the optically- 
thick limit (with E* (x 'S,g/K) is currently difficult to test and 
requires further high resolution observations. It is important 
to stress that in our models the star formation efficiency in the 
optically-thick nuclei of starbursts is actually lower than what 
would be inferred from an extrapolation of the star formation 
efficiency in normal star-forming galaxies because of the extra 
support from radiation pressure in the optically-thick limit. 

Our calculations show that in the optically-thick portion of 
starburst disks near ^100 pc, the flux, star formation rate 
per unit area, and the effective temperature of the disk are 
roughly constant with characteristic values of E ~ 10^^ L© 
kpc“^, S* ~ 10^ Mq yr“^ kpc“^, and T^ff ~ 90 K (§2.3). To test 
these predictions, we have estimated the fluxes (temperatures) 
in the emission regions of ULIRGs using the resolved radio im¬ 
ages of Condon et al. (1991) as indicative of the sizes of nuclear 
starbursts in these systems (the fluxes inferred using FIR black- 
body temperatures may not accurately characterize the fluxes 
in the emission region if there is still significant obscuration at 
>20-30 microns). We find excellent agreement with the pre¬ 
dictions of our models (Figs. 3 & 4). This supports our interpre¬ 
tation that a significant fraction of the radiation from ULIRGs 


is produced by an Eddington-limited starburst. In the future, 
imaging at optically-thin FIR wavelengths will test our assump¬ 
tion that the radio emission traces star formation in ULIRGs. 

In our model of starburst disks, the central temperature ex¬ 
ceeds the dust sublimation temperature (Tjub ~ 1000 K) at radii 
~ 1 - 10 pc. The opacity of the disk then drops dramatically 
and the star formation rate required to maintain Q ^ 1 is very 
large (Fig. 2). This poses a severe problem for fueling a central 
AGN with gas stored in the starburst disk at larger radii (Sirko 
& Goodman 2003). At first glance it appears that all of the 
gas is converted into stars in the “opacity gap” with little left 
available to power a central AGN. 

A number of solutions to this problem have been proposed, 
all of which hinge on the efficiency of angular momentum 
transport in the disk (Shlosman & Begelman 1989; Goodman 
2003). For fixed M, higher viscosity implies lower surface den¬ 
sity and thus self-gravity is comparatively less problematic. For 
this reason, we have focused on the possibility that angular mo¬ 
mentum transport proceeds via a global torque, as would be 
provided by a bar, spiral waves, or a large-scale wind. With this 
prescription for angular momentum transport we have extended 
our starburst disk models to smaller radii where they connect 
consistently with AGN disks on sub-parsec scales. A key in¬ 
gredient to these calculations is an equation that accounts for 
the radial loss of gas due to star formation (eq. [43]). This 
limits the star formation rate required to support the disk in 
the opacity gap region and, therefore, allows for significant gas 
accretion and AGN fueling (Appendix A). In this respect our 
model differs from that of Sirko & Goodman (2003), who also 
considered star formation-supported AGN disks, but assumed 
constant M. 

Our model provides a way of quantifying the connection be¬ 
tween starbursts and AGN activity. In particular, we find two 
classes of disk models (Figs. 5-7): (1) For mass supply rates at 
the outer edge of the disk less than a critical rate Me (eq. [44]), 
a starburst occurs primarily in a narrow ring at Rout- Nearly all 
of the gas is converted into stars at large radii and the spectrum 
is dominated by a FIR starburst peak. This class of solutions 
corresponds to the limit in which the star formation time in the 
disk is shorter than the viscous time (r* < Tadv)- (2) On the 
other hand, forMout > Me, a fraction of the gas accretes inwards 
to smaller radii because the advection timescale is shorter than 
the star formation timescale. These disk solutions have an outer 
starburst at ^ Rout and an inner nuclear starburst ring on ~ 1 -10 
pc scales at the opacity gap (Fig. 5). In addition, for fiducial pa¬ 
rameters we find that the accretion rate onto the central BH is 
~ 1 - 10 M© yr“^ (Appendix A). Models with Mout > Me can 
be dominated by AGN emission, although there is a prominent 
contribution in the mid- and far-IR from the inner and outer 
starbursts (Fig. 7). Both the broad but sub-dominant IR peak 
and the magnitude of the AGN luminosity in our models are 
reasonably consistent with composite AGN spectra. 

One prediction of our models is that a significant fraction of 
the IR emission in luminous AGN must be from star formation 
rather than reprocessing. The star formation is required to sup¬ 
port the accretion disk at large radii against its own self-gravity 
(see Rowan-Robinson 2000 for an observational discussion of 
this point). We have identified an additional observational con¬ 
sequence of the parsec-scale starburst that is predicted to occur 
when the temperature of the disk reaches the sublimation tem¬ 
perature of dust. Our estimates suggest that for luminous AGN, 
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the nuclear starburst is able to inflate the photosphere of the disk 
to a height hpi, ~ r (§3.3 and Fig. 8). This extended atmosphere 
of the disk contains little mass but produces significant obscu¬ 
ration, and may account for some of the nuclear obscuration 
observed in Type 2 AGN. 

As a final application of our calculations, we note that there is 
direct observational evidence for compact starbursts near mas¬ 
sive black holes. In particular, observations of the Galactic Cen¬ 
ter reveal a population of young O and B stars within the central 
parsec of our Galaxy. Most of these stars appear to lie in a thin 
disk located within 0.1 -0.3 parsecs of the black hole (Levin 
& Beloborodov 2003; Genzel et al. 2003). The origin of these 
stars remains uncertain, but their disk-like kinematics suggests 
that they might have formed in a dense self-gravitating accre¬ 
tion disk, if the black hole in the Galactic Center were accreting 
at a much higher rate several million years ago (e.g.. Levin & 
Beloborodov 2003). To explore this hypothesis in the context 
of our models. Figure 9 shows the star formation rate and accre¬ 
tion rate as a function of radius in models with a = 75 km s~*, 
Mbh = 4 X IO^Mq, /?out = 3 pc, and an a-viscosity with a = 0.3 
(we also increased the opacity by a factor of three relative to 
the curves presented in Figure 1 to account for the super-solar 
metallicity of the Galactic Center). The chosen outer radius is 
appropriate if the circum-nuclear disk in the Galactic Center — 
a current reservoir of ~ 10“*- 10^ Mq of gas (e.g., Jackson et 
al. 1993; Shukla, Yun, & Scoville 2004) — were accreted onto 
the central black hole. The mass supply rates we have used 
— Moui = 0.015 - 0.15 Mq yr“^ — are in the range expected if 
the circum-nuclear disk were accreted on a viscous timescale 
~ 10® years. Figure 9 shows that the star formation rate in this 
hypothesized burst of accretion typically has a strong peak at 
r ~ 0.1 pc. This is the “inner nuclear starburst” caused by the 
decrease in opacity when the central disk temperature exceeds 
the sublimation temperature of dust. The location of this star- 
burst is strikingly similar to the current location of the stellar 
disk in the Galactic Center.^® 

Several aspects of our calculations require further study. The 
first is the stability of our disk solutions. It is well known that 
radiation-pressure supported disks are prone to numerous insta¬ 
bilities: viscous, thermal, convective, and photon-bubble (e.g., 
Lightman & Eardley 1974; Piran 1978; Blaes & Socrates 2003). 
Such instabilities have been extensively studied in the context 
of the central ~ 10-100 Schwarzschild radii of black hole ac¬ 
cretion disks, where the radiation pressure is on free electrons. 
An important difference between our solutions and these is that 
at large radii the disk is heated primarily by star formation, not 
accretion. We find that as a result, the disk is globally ther¬ 
mally and viscously stable (Appendix B). As in the local ISM, 
however, the disk may fragment into a multi-phase medium. A 
second aspect of our model that requires further study is the 
physics of star formation under the unusual conditions appro¬ 
priate to gravitationally unstable accretion disks, particularly 
at radii ~ pc where the gas is much denser and hotter than in 
normal star-forming regions. Because Q <C 1 in the absence of 


star formation, and because the cooling time of the disk is much 
shorter than the orbital time (Appendix B), we believe that frag¬ 
mentation of the disk is inevitable. But beyond this, the prop¬ 
erties of star formation under these conditions are poorly un¬ 
derstood (Goodman & Tan 2004 argue that the disk will form 
supermassive stars). Interestingly, the power-law solutions for 
the star formation rate in the opacity gap region (Appendix A) 
are independent of the IMF-dependent parameter e (eq. 10). On 
larger scales (r > 1-10 pc), however, a top-heavy IMF (larger 
e) can suppress the star formation rate (eqs. 16 & 28), decrease 
the star formation efficiency rj (eqs. 18 & 30), and decrease the 
critical mass accretion rate required for AGN fueling (eq. 44). 
Another component of our model requiring further study is the 
interaction of the AGN emission with the flared starburst disk 
and the tenuous dusty obscuring atmosphere predicted in §3.3. 

Finally, we note that large-scale outflows of cold dusty gas 
are commonly observed from starbursting galaxies, including 
ULIRGs (Alton et al. 1999; Heckman et al. 2000; Martin 2004). 
We have recently proposed that radiation pressure on dust is an 
important mechanism in driving such winds (Murray, Quataert, 
& Thompson 2005 [MQT]). The disk models presented here 
provide support for this idea because it is natural to suspect 
that a radiation pressure supported starburst disk will drive a 
wind from its photosphere, much as near-Eddington stars drive 
strong outflows. In MQT we further argued that both star- 
bursts and AGN have a maximum Eddington-like luminosity 
given by Lm = 4-fga‘^c/G}^ One can show that the disk mod¬ 
els presented in this paper have luminosities < Lm, with the 
equality obtained (approximately) when /r ~ r and the spheri¬ 
cal limit employed in MQT is realized. We also note that the 
results in this paper, together with those of MQT, imply that 
fueling a luminous AGN requires an accretion rate in the range 
Me < Mout ^ 47niax = With the black hole accretion 

rate given by equation (A4), this yields a ratio of star formation 
rate to black hole accretion rate of ~ 100- 1000, similar to the 
observed ratio of stellar mass to BH mass in nearby galaxies 
(Magorrian et al. 1998; Haring & Rix 2004). 
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For the models in Figure 9 that produce a peak in M* at r ~ 0.1 pc, the continuum approximation (see footnote 13) breaks down for r < 0.05 pc. In this region the 
disk will not self-regulate and it will cool rapidly and form more stars. The star formation rate in this region will then be determined by the criterion r//r ~ Nm* and 
not 2 ~ 1- Because the continuum limit is violated by a factor of ~ 3 —5 in the models displayed in Figure 9 for r < 0.05 pc, M* will be 3 —5 times higher in this 
region. The gas accretion rate at yet smaller radii will be only somewhat smaller. 

The “M” in Lm stands for “momentum-driven”. 
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APPENDIX 


APPENDIX A: THE BLACK HOLE ACCRETION RATE 


An interesting feature of the results shown in Figure 5 is that all of the models with Mout > Me deliver precisely the same amount 
of mass to the AGN disk at small radii. More specifically, all of the dashed lines for Mout > 220 Mq yr“' in Figure 5 come together 
and attach to a single power-law solution in the opacity gap between r ~ 1-10 pc, and then settle onto a constant M « 4 Mq yr“^ 
solution at smaller radii. In this appendix we explain why this is the case. 


The steep drop in the opacity in Figure 1 occurs when the central disk temperature reaches the sublimation temperature Tjub at a 
radius /?sub- In the models presented in Figure 5 with M > Me, Rsub corresponds to the sudden increase in M* between ~ 1 and ~ 10 
pc. The rapid decrease in k inside the opacity gap can be approximated by a power law, k(T > Tsub) = n&ubT^■ For example, in the 
opacity models of Bell & Lin (1994), /3 = -24 and T^^b ~ 900 K (see Fig. 1). Because the star formation rate required to maintain 
2 ^ 1 in the optically-thick limit is given by S* ex (eq. [28]), one might expect E* to increase by ~ 5 orders of magnitude 
inside the opacity gap because k decreases by this amount at T^^b (see Fig. 1). In the models with constant fg shown in Figure 2, 
this is precisely what happens. However, when we self-consistently solve for M(r) (eq. [43]), as in the models presented in Figure 
5, the results are qualitatively different. The star formation rate increases by only ~ 1 order of magnitude at r ~ R^^b because of an 
important physical effect. As k decreases and M* increases, the accretion rate must decrease because the gas is depleted as a result of 
star formation. This leads to a lower surface density and, therefore, a lower star formation rate is required to maintain Q ~ 1. Thus, 
an equilibrium is established between star formation and advection in the opacity gap (radv = t*). An analytic solution representing 
this balance between star formation and accretion can be derived by solving equation (43) for M(r) using k cx T^, T (X (eq. 
[26]), S* cx Eg/K (eq. [28]), and the fact that there is a one-to-one relationship between fg and M. If angular momentum transport 
proceeds via global torques, the fg-M relation is given by equation (46), and the analytic solution for M(r) takes the form 

M(r)'/2+^/4-M(/?,ub)‘/2'"^/4 =A (^r^/4+5/3/8_^3/4+5/3/8^ ^ 

where A is a /3-dependent constant that will be specified below. Because /3 is <C 0 in the opacity gap, the solution quickly loses 
memory of the initial conditions at R^ub- As a result, a single universal equilibrium power law solution for M « M* exists in the 
opacity gap. Because of this equilibrium between star formation and advection, the increase in M* at /?sub is much less pronounced 
in these solutions than in the constant fg models (compare Fig. 2 & 5). 


The general expression for the constant A for arbitrary /3 is a bit awkward, but the limit (3 —> -oo (as appropriate for a step function 
in opacity) provides a useful approximation. In this limit, the accretion rate inside the opacity gap is given by 
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where Fwoo = Tsub/IOOOK, mo.i = m/0.1, rw = r/lOpc, and Mg = Mbh/IO^Mq. Note that the normalization of M in equation (A2) 
depends only on the black hole mass Mbh, the Mach number m = Vrjcs, and the dust opacity law through T^^b- This explains why all 
of the solutions with Mout > Me in Figure 5 are the same inside the opacity gap. Note also that because M(r) is a power-law inside 
the opacity gap, the star formation rate M*(r) is as well (see eq. [43]). Indeed, for these solutions, M*(r) = Trr^E* = (5/4)M(r), in 
good agreement with Figure 5. Finally, because M* cx in the opacity gap, both 7/^ and E* are proportional to The central 
temperature T and the surface density Eg, on the other hand, are independent of radius. 


As is clear from Figure 5, the equilibrium power law solution for M given in equation (A2) does not persist to arbitrarily small 
radii. Because p increases and T is constant in the opacity gap, gas pressure eventually dominates radiation pressure. The two 
pressures are comparable at a radius /?eos given by 

Reos ~ 1 MfMfl^niilQr^l‘^ pc, (A3) 

where M\ = M/1 Mq yr“'. For r < /?eos, one can show that star formation is no longer a significant impediment to accretion: M* 
drops precipitously and M approaches a constant. As a result, we can estimate the accretion rate onto the black hole by evaluating 
equation (A2) at r ~ 7?eos- This yields 

Mbh = 2^/'7;yymMyyG-i/2 (^) j yr'*, (A4) 

where we emphasize that m is the Mach number in the opacity gap region; even if angular momentum transport becomes inefficient 
for r < Reos (e-g- m <C 0.1 or a transition to local viscosity with a < 0.1 occurs), Mbh remains unchanged. Using equation (A4) and 
taking Lbh = CMbhc^ we find that 

Lbh ~ lO^^Co.iergs s'^, (A5) 


where ^o.i = C/O L This result is in excellent agreement with the models presented in Figures 5 and 7. With Ledd = (Meddc^ = 
4ttGMbhc/ Kes, where Kes is the electron scattering opacity, we find that Mbh/Medd ~ O.lmo.iCo.iMg^^^. Therefore, for typical 
parameters, the black hole is fed at a reasonable fraction of its Eddington rate. Because of the scaling with Mbh we expect black 
holes of smaller mass to be preferentially super-Eddington. 
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Note again that from equation (A4), the accretion rate through the opacity gap into the AGN disk depends primarily on the dust 
sublimation temperature and on the rate of angular momentum transport in the opacity gap (via the Mach number m). Because of 
the latter dependence, it is of interest to compare the above expressions with the analogous results assuming that angular momentum 
transport in the opacity gap is via a local viscosity, for which the fg-M relation is given by equation (42). In this case the accretion 
rate in the opacity gap is given by 

M(r) = (^^7;4,V27rGe) ' ^ (A6) 

and the accretion rate onto the black hole is given by 

Mbh = ^ « 2 X lO-3ao.i0“'riYooMe yr"'. (A7) 

GQ \ nip / 

Equations (A6) and (A7) highlight the fact that in our models, local viscosity is insufficient to transport matter through the opacity 
gap at a rate capable of fueling a bright AGN. 


APPENDIX B; disk STABILITY 

Appendix C lists the full set equations we solve for calculating the disk models presented in this paper. Because of the complicated 
temperature dependence of the opacity (Fig. 1), together with equations (C3) and (C5), we find multiple solutions at some radii. These 
solutions have the same pressure (by eq. [C3]), but different temperatures and opacities. In order to determine whether all of these 
solutions are physical, we assess the thermal and viscous stability of our disk solutions in this Appendix. A second motivation for 
doing so is that radiation pressure dominated disks close to black holes are known to be thermally and viscously unstable in certain 
circumstances (e.g., Piran 1978), and so it is important to check whether the same is true for our solutions. 


B.l Timescales 


There are four relevant timescales characterizing the disk: (1) the dynamical timescale, Tdyn ~ (2) the advection timescale, 

Tadv = r/Vr, (3) the star formation timescale, r* = (pH)"', and (4) the thermal or cooling timescale, Tth = J^gC^/F, where F is the flux. 
The latter is also the photon diffusion time across the disk. 


In the optically-thick limit, and assuming that radiation pressure dominates gas pressure (prad/Pgas ^1), the ratio of the cooling 
timescale to the dynamical timescale is very small: 


Tth 

"^dyn 
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For the disk solutions presented in Figure 5, the optical depth is always greater than unity, but not larger than ~ 1000, even for r ~ 1 
pc. Typically, the optical depth is closer to ~ 10- 100. This, together with the fact that q ~ fgir (eq. [6]) implies that Tth/rdyn 1. 
This is an important difference relative to canonical accretion disk models. For a local a viscosity an accretion-heated disk has 

Tth ^ tr Tdyn ^ Tdyn. 


We may also compare the star formation timescale with the dynamical timescale. Again assuming that Prad/Pgas ^ 1 and that 
Tv > 1, we find that 


^dyn 



(B2) 


As a consequence, throughout our disk models T*/Tdyn > 1. This result is, of course, in keeping with our assumption that star 
formation occurs on a timescale longer than the local free-fall timescale (that is, p < 1). 


Fastly, we compare the advection timescale with the dynamical timescale. If angular momentum transport is produced by a global 
torque (as in Figure 5), Tadv ~ Tdyn(T//i)m“', which is always greater than unity. If angular momentum transport is instead driven by 
local viscosity, (Tadv/Tdyn) ~ (r/lif’aT^. As discussed in §3, the advection timescale can be larger or smaller than the star formation 
timescale, depending on the ratio of M to Me- 


These considerations lead to the following well-defined hierarchy of timescales in Q ^ 1 disks: 

Tth < Tdyn < T*. (B3) 


B.2 Thermal Stability 

The fact that Tth is much less than Tdyn has important implications for the thermal stability of the disk. In particular, because of 
this disparity of timescales, the star formation rate per unit area E* is constant on a timescale Tth. This implies that the volumetric 
heating rate q'^ ^ F /hoc. E*//i is a constant on Tth, and, hence, d\nq^/dlnT = 0. However, the volumetric cooling rate in the optically 
thick limit is q~ ~ /(jy/h). Taking k oc T^, d\nq'^/dlnT oc (4-/7). Therefore, if /3 > 4, the solution is thermally unstable. 

This condition is only obtained when k cx at T ~ 2000 K (see Figure 1; Bell & Fin 1994). All other solutions are thermally 
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stable. In particular, our 2 ~ 1 radiation-pressure supported starburst disk with T ~ 100s K is globally thermally stable, in contrast 
to radiation-pressure supported disks heated by accretion energy. 

In the optically thin limit, the cooling rate is q~ ^ j ^ kB (x where j is the emissivity. Therefore, if /3 < -4, the solution is 
unstable. There are only two regions where such a condition might obtain and they correspond to sudden decreases in opacity (nearly 
step functions in k[T]) at T ~ 100-200 K and in the opacity gap at T ~ 1000 K in Figure 1. Otherwise, all optically-thin solutions 
are globally thermally stable. 


B.3 Viscous Stability 

To evaluate the viscous stability of our disks we must compute the partial derivative dM/d'Sg (e.g., Lightman & Eardley 1974). 
For dM/d'Sg < 0, the disk is viscously unstable and we expect perturbations to grow and disrupt the disk on a timescale Tadv- If 
angular momentum transport is driven by global torques, M (x EgCj. But as a consequence of our assumption of marginal stability 
against self-gravity (Q ~ 1), at any radius, Cj oc oc Moreover, because Tadv ^ Tdyn it is reasonable to maintain Q ~ 1 during 
perturbations to the surface density. These results imply that dM/d'Sg is always greater than zero and the disk is viscously stable. A 
similar argument applies if angular momentum transport is instead driven by local viscosity (M oc ~ 1 disks thus appear to 

be viscously stable. 


B.4 Physical Disk Solutions 

In the region where Q ~ 1, we often find three solutions to our disk equations, one low-temperature optically-thick solution and 
two high-temperature optically-thin solutions (see also Sirko & Goodman 2003). Throughout this paper we have focused on the 
low-temperature optically-thick solution as the physical solution. Here we justify this choice. 

All three of the of the disk solutions appear viscously stable. The two high-temperature solutions bracket the opacity gap: one has 
T = Tsuh ~ 1000 K and the other occurs at T ~ 2000-5000 K. By the optically-thin thermal stability criteria of §B.2, the solution 
with T ~ 1000 K is thermally unstable because f3 is very large and negative. This solution is therefore probably unphysical. The 
solution with T ~ 2000-5000 K is formally thermally stable because f3 is positive (Fig. 1). However, the opacity curve employed 
here does not take into account important line cooling processes in this temperature range. Indeed, studies of the ISM have shown 
that gas with T ~ 2000-5000 K is thermally unstable (Wolfire et al. 1995). Although the ISM of dense starburst nuclei may be 
interestingly different than that of normal star-forming galaxies, it is still likely that optically thin gas at several thousand Kelvin is 
thermally unstable. For this reason we do not believe that the T ~ 2000-5000 K solution represents a physical global equilibrium of 
the disk. It is, however, likely that the existence of multiple solutions at the same pressure implies that the ISM is prone to breaking 
up into a multi-phase medium. A detailed investigation of the multi-phase structure of the ISM in starburst galaxies is clearly of 
interest, but beyond the scope of this paper. 


APPENDIX C: DISK EQUATIONS 


The equations employed in constructing the models shown in Figures 5-7 and described in §3 are 


H(r) = HK(r) = 


/ GMbh ^ 
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In the outer part of the disk where accretion heating is insufficient to maintain 2 ~ 1, these equations are solved subject to the 
conditions 

P = 


V2ttGQ 


(CIO) 
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(with Q= 1) and 

asBT^ff=yAc^+^M[l-V^)n\ (Cll) 

In the inner part of the disk where accretion heating maintains Q> equations (CIO) and (Cl 1) are replaced by 

(l - v/^) (C12) 

The radius at which Q > 1 is denoted /?agn- Throughout this work we set ^ = 1 in equation (C3) and e = 10“^ in equation (Cll). 
Mbh is specified by Mbh = 2 x lO^u^oo The inner edge of the AGN disk is taken to be at = 3(2GMbh)/c^. A7out and /?out are 
the mass supply rate at the outer boundary and the location of the outer boundary, respectively. In equation (C8), we generally take 
m = 0.1 —0.2. 


When calculating disks whose accretion is driven by a local viscosity, as in Figure 9, we replace equation (C8) with 


M = liivYig 


c/lnfl 


c/lnfl 

d\nr 
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(Cl 3) 


APPENDIX D: the critical MASS SUPPLY RATE 

As shown in Figure 5 and Appendix A, when Mout exceeds a critical rate Me, gas can accrete to small radii fueling a bright AGN. 
By contrast, for < Me, star formation consumes the gas in a narrow range of radii at ^ /?out- In the optically-thick limit, and to 
the extent that n may be approximated by kqT^, the critical mass supply rate is given in equation (44). 

The balance between star formation and accretion is different in the optically-thin limit than it is in the optically-thick limit. To 
show this, we note that equation (43) admits a simple power-law solution when angular momentum transport is driven by a global 
torque (eq. [45]): 

M(r) (X (Dl) 

Because 

-7(T3oo(e3Cmo,i)“' (D2) 

V2e^cm 

we see that unless there is a very efficient global torque (m > 1), gas will be converted into stars rather than accreting to smaller radii. 
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Fig. 1. — The Rosseland mean opacity as a function of temperature at several densities using two opacity models from the literature. Note the pronounced decrease 
in the opacity at the sublimation temperature of dust at T 1000 K (the “opacity gap”). The rise in the opacity at higher temperature is due first to H-scattering, 
then to bound-free and free-free interactions, and finally to scattering off of free electrons. The thick solid line shows that at low temperatures, the opacity scales 
with the square of the temperature. 
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Fig. 2.— Vertical optical depth (log^Q[7v]; upper left hand panel), opacity (logjQ[/^ (cm^ g~^)]; lower left hand panel), star formation efficiency (log|Q[? 7 ]; lower 
right hand panel), and both the star formation rate per unit area and the effective temperature (logjQ[S* (MQyr“^ kpc“^)], logjQ[reff (K)]; upper right hand panel) 
as a function of radius for models with u - 200km s“^ and either fg - 0.03 (dashed line) or /^ = 1 (solid line), r] approaches a constant for r > 100 pc in the 
model with fg = 0.03 when the disk is optically thin (see also eq. [18]). The model with fg = \ demonstrates that the star formation rate per unit area is roughly 
constant at intermediate radii 100 pc (see §2.3). At small radii ~ 0.01 —10 pc, S* is extremely large in the region where the dust sublimates and the opacity drops 
precipitously. Such large star formation rates are unphysical since the mass accreting through the starburst disk would quickly be exhausted. In §3 we construct 
more realistic models, taking into account the depletion of the gas locally as a result of star formation. Including this effect significantly reduces S* on few-parsec 
scales (compare with Fig. 5). 
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Fig. 3.— Histogram of the number of ULIRGs at a given flux using the sample in Tables 1 & 2 of Condon et al. (1991); the bins have a width of 0.2 in log^Qp. 
We estimate the intrinsic flux in the nuclear starburst using F = LpiR/Ai-adio^ where Aj-adio is the area of Condon et al’s elliptical Gaussian fit to the radio emission; 
we exclude the 4 unresolved sources. The peak in the histogram at 10^^ Lq kpc“^ is in good agreement with the predictions of our optically thick disk models 
(§2.3). See the text in §2.5 for uncertainties associated with using the radio size as a proxy for the starburst size. 
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Fig. 4.— Flux as a function of radial scale for the local ULIRGs in Condon et al. (1991) (filled squares). The fiux is inferred using F = TpiR/Aradio^ where Aradio 
is the area of Condon et al’s elliptical Gaussian fit to the radio emission; we define the radial scale R by Aradio = The solid lines show results from our constant 
fg models (compare with Figs. 2). From lowest to highest flux the models have a = 200 km s“^ and fg = 0.1, 0.3, and 1.0 (solid lines). For comparison, the dashed 
line shows a model computed with a = 300 km s“^ and fg = l .0. 
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Fig. 5.— The local star formation rate (M* = solid lines) and accretion rate (dashed lines) as a function of radius for Mout = 80, 160, 220, 320, and 640 

Mq yr“^ in a model with Roui = 200 pc and a = 300 km s“^. Note the strong bifurcation between models with low mass supply rates (Mout < 160 Mq yr“^; see 
eq. [44]) and models with high mass supply rates (Mout ^ 220 Mq yr“^). The latter produce outer (r Rout) and inner (r ~ 1-10 pc) starbursts, but also fuel a 
bright central AGN with Mbh ~ 4 Mq yr“^ (Appendix A). Figure 7 shows the computed spectra for each model. Figure 6 shows the full disk structure for the 
model with A/out = 320Mq yr“^ . 






logio[R(pc)] log,o[R(pc)] log,o[R(pc)] 


Fig. 6. — logic of T (solid line, upper left panel), Tgff (dashed line, upper left panel), p (upper middle panel), (upper right panel), radiation pressure pr (solid 
line, middle left panel), gas pressure pg (dashed line, middle left panel), f3 = pg/ptot (middle panel), h/r (middle right panel), Cs (lower left panel), k. (lower middle 
panel), and rv (lower right panel) all in cgs units as a function of radius in the model with Mout = 320 Mq yr“^, Rout = 200 pc, and m = 0.2. The opacity gap region, 
highlighted in Figure 5, is located at r 1 pc. 
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Fig. 7.— The spectral energy distribution as a function of wavelength for each of the disk solutions presented in Figure 5. Solutions with small mass supply rates 
(Mout < 160 Mq yr“^ < Me) are dominated by a starburst at large radii that produces a prominent FIR peak. Solutions with large mass supply rates (Mout > 220 
Mq yr“^ > Me) have a similar FIR bump, but also have a spectral peak at ~ 10 ^m due to the inner ring of star formation in the opacity gap at 1 -10 pc. As shown 
in Figure 5, these solutions also provide mass to the black hole at a rate comparable to the Eddington limit. The total bolometric emission from these models can be 
dominated by the central AGN, although there is a significant contribution from star formation (Appendix A). 
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Fig. 8 .— Estimates of the midplane (solid lines) and photospheric (dashed lines) scale-heights for two models with Mout = 320 and 640 Mq yr ^ The model 
with Mout = 640 M 0 yr~^ has higher midplane hjr aX Rout and higher photospheric hp^jr aX r 10 pc. Both models have the same midplane and photospheric scale 
height for r < 2 pc because of the power law solution in the opacity gap (Appendix A). Note that for 0.1 0.5 pc the disk is gas pressure dominated at the 

midplane, but radiation pressure dominated at the photosphere (see Fig. 6 ). Values of hphjr > 1 indicate that the photosphere is unbound. 
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Fig. 9.— Star formation rate (solid lines) and gas accretion rate (dashed lines) in disk models appropriate to the Galactic Center (Mbh = 4 x IO^Mq and 
u = 75 km s“^). We consider a mass supply rate Mout at a radius Roni = 3 pc and accretion via a local a viscosity with a. = 0.3. The chosen values of i?out and Mout 
are motivated by observations of the circumnuclear disk in the Galactic Center. Note that the opacity was increased by a factor of 3 relative to Figure 1 to account 
for the super-solar metallicity in the Galactic Center. These results show a ring of star formation at 0.1 pc where there is currently a young disk of O and B stars. 




